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Abstract 



We present an analytic and numerical analysis of the Gledzer- 
Ohkitani-Yamada (GOY) cascade model for turbulence. We concen- 
trate on the dynamic correlations, and demonstrate both numerically 
and analytically, using resummed perturbation theory, that the cor- 
relations do not follow a dynamic scaling ansatz. The basic reason 
for this is the existence of a second quadratic invariant, in addition to 
energy. This implies the breakdown of the Kolmogorov type scaling 
law, in a manner different from the conventional mechanisms proposed 
for Navier- Stokes intermittency. By modifying the model equation so 
as to eliminate the spurious invariant, we recover to good accuracy, 
both dynamic scaling and the Kolmogorov exponents. We conclude 
that intermittency in the GOY model may be attributed to the effects 
of the spurious invariant which does not exist in the 3-dimensional 
Navier-Stokes flow. 
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1 Introduction 



A class of models of liydrodynamic turbulence that has recently attracted 
some interest are the so called 'shell' or Gledzer-Ohkitani-Yamada (GOY) 
[1][2] models. These models exhibit numerically a clear deviation from the 
expected scaling predicted by an application of Kolmogorov like dimensional 
analysis. In addition, high moments of the variables exhibit muUiscaling. 
It has been speculated that the mechanism for generation of multiscaling in 
these models may be related to the possibility of multiscaling in Navier-Stokes 
turbulence. Intermittency in hydrodynamic turbulence is usually attributed 
to large spatial fluctuations in the local energy dissipation e. In the shell 
models there is no space coordinate, and therefore it is unclear how this 
mechanism could be applied to shell models. However, reference [2] considers 
the connection between dissipation fluctuations in time and intermittency in 
the structure function via a 'Taylor hypothesis' type argument. 

We demonstrate in this paper that the "intermittency" observed in the 
GOY models is largely and perhaps completely due the existence of a second 
quadratic invariant in addition to energy [3] [4]. This quadratic invariant 
does not have an analog in Navier-Stokes equations, and is not even positive 
definite, but still its existence will be shown to have important effects on 
the dynamics. This claim will be supported by comparison with a modified 
version of the GOY model that lacks the second invariant. In particular the 
original GOY does not possess dynamic scaling, while the modified version 
obeys dynamic scaling. 

The paper is organized as follows; Section 2 defines the model and de- 
tails some of its important symmetries and conserved quantities. In section 
3 we present our numerical results and demonstrate the absence of dynamic 
scaling in the GOY model. In section 4 we examine the direct interaction 
approximation (DIA) equations for the GOY model and conclude that dy- 
namic scaling is indeed inconsistent with the DIA due to the existence of 
the second invariant. In section 5 we consider the modified GOY and show 
that dynamic scaling and the Kolmogorov exponents are restored. Section 
6 offers some conclusions regarding the relations between intermittency in 
shell models and Navier-Stokes turbulence. 

2 General properties of the model 
2.1 Model definition 

The GOY model reads (in its general form) [I], 

[dt + V kl)Un = l[a kn Un+lUn+2 + b k^-l Un-lUn+l+C kn-2 Un-lUn-2]* + fn, (1) 

where n > is the complex dynamical variable of shell n, with wave 
number kn = koq^. q > 1 is the shell spacing, p is the 'viscosity', is the 
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external forcing of shell n, and a, b and c are real coupling constants whose 
sum is zero. The number of shells may in principle extend to infinity but in 
practice the presence of the viscosity term acts as a cut-off. This allows us 
to truncate the system of equations after a finite number of levels without 
changing the dynamics. 



2.2 A phase symmetry 

The unforced Navier Stokes equation possesses translational invariance. In 
k space, this implies an invariance of the equation with respect to the trans- 
formation Vk — )• e^^'^°V]^ of a solution . This symmetry implies that the 
Fourier components of the velocity field are delta correlated in A;, e.g. 

(vkVq) = F(k),5(k + q). (2) 

This follows from the assumption that the symmetry of the equation is not 
broken by the forcing. The implication of this relationship is that the pair 
correlators and the Green's functions are all diagonal in A;, in other words, 
the vertex conserves momentum. 

The translational invariance of the Navier-Stokes is not retained in the 
GOY model. However, there exists an analog of this symmetry which leaves 
the GOY model equations invariant. This is the transformation 

Un e^'^Un (3) 

Un+1 ^ e'^^-'^u^+^ (4) 

Un+2 e~"^Un+2 (5) 

Note that 6 and (f> are arbitrary global phases. This transformation adds the 
same phase to shell variables that are separated by a distance of three shells. 
A reasoning similar to that of the Navier-Stokes case leads to the conclusion 
that 

= if m is not divisible by 3. (6) 

This symmetry also implies that the shell variables have a zero mean. Similar 
relations may be obtained for higher order correlations. 



2.3 Quadratic conserved quantities 

A crucial ingredient of the Kolmogorov argument is the existence of one 
(quadratic) invariant for the Euler equations, energy. For the Navier-Stokes 
equations The energy is dissipated at the smallest scales and is injected at 
the largest scales. Kolmogorov assumed the existence of an energy fiux that 
is scale independent in the inertial range. Motivated by this, it is useful to 
analyze the GOY model in terms of quadratic invariants. 
We are interested in conserved quantities of the form 

C = J2A^\u^\' (7) 
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and demand that in the hmit of no forcing or viscosity dfC = 0, independent 
of Un- Hence, inserting in equation (1) we obtain a system of equations 

Ana + An+lb+ An+2C=0. (8) 

This system may be viewed as hnear recursion relation for An+2 in terms of 
An- This means that there are only two linearly independent solutions 
for the A's. As we demand energy conservation, one of these solutions is 
An = 1 and consequently a + 6 + c = 0. We look for the second independent 
solution in the form An = A". We may solve for A yielding A = a/c. 
With the standard choice of parameters a = l,6 = c= — 1/2 we obtain 
An = (—2)" = [ — )^kn. This new conserved quantity L = is 
spurious in the sense that it has a dependence on the parameters of the model. 
In the usual choice of parameters (see eq (12) below) this invariant is not 
positive definite, however this does not prevent it from having an important 
dynamical effect [3] [4] [5] (see below). For a particular set of parameters, it is 
possible to interpret L as the analog of helicity which is conserved in the Euler 
equations [3] . However the existence of this invariant was not considered as 
a requirement in the original definition of the model [1]. 

3 Numerical examination of scaling proper- 
ties 

In previous works [2] [6] it has been found that the simultaneous correlators 
of the shell variables obey scaling laws, such as 

(i^.n ~ (9) 

for shell numbers inside the inertial range, (q was found to have a non-linear 
dependence on q. Such a behaviour is termed 'multiscaling'. The inertial 
range is defined as the shells in which both the dissipative and the forcing 
terms are small with respect to the non-linear term. 

The usual Kolmogorov picture implies also that the different-time corre- 
lators should obey a dynamic scaling relation of the type 

{u:it)unit + T))r.ki^f(rK). (10) 

/ is a universal scaling function independent of the specific shell. This scaling 
relation becomes in the frequency domain 

{\Uniu)\')r.kyjiu/K), (11) 

where y + z = (2. Kolmogorov scaling predicts that y = —4/3 and z = 2/3. 
Although intermittency may change the values of these exponents, one would 
still expect naively the dynamic scaling hypothesis to hold. 
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In order to test this hypothesis we integrated numerically the GOY equa- 
tions and considered the frequency power spectra of the shell variables. Fol- 
lowing the usual procedure we chose /s = (1 X 10""^ and all other /„'s to 
be zero. In the numerical integration of the equations we used the following 
parameters: 

a = l, 6 = -0.5 c=-0.5 , , 

r- r- 12 

A;o = 1/(32^2) u = 6 x 10'^ q = ^2. 

There were 43 shells (running from to 42) in the model which corresponds 
to roughly 21 level in the q = 2 models. The equation were integrated using 
a 4th-5th order Runge-Kutta adaptive time step method. The integration 
time was 6.6 X 10"* which corresponds to approximately 1500 forcing scale 
turnover times. 

The power spectra were calculated as follows: A data string of 65536 
equally time-spaced samples, calculated by linear interpolation of the values 
obtained from the integrator, were fast Fourier transformed with a Hanning 
window, and |m„(c<;)P was calculated. This procedure was repeated many 
times and the values of |m„(c<;)P were averaged. 

The spacing between the samples was chosen to be somewhat smaller 
than the characteristic time scale for the highest shell (the 37th) which was 
analyzed. The size of the sample was sufficient to capture the whole dynamic 
range of the power spectrum except for the lowest shells. 

An attempt to collapse the data of several power spectra is shown in 
figure 1. The power spectra are rescaled according to the typical height and 
width of each power spectrum as follows: 

fi^\w) = ^U{wu^), (13) 



where 
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w is the dimensionless frequency io/iOn- 

The three power spectra shown in figure 1 demonstrate that the rescaling 
performed in (14) yields a poor data collapse. There is a general trend for the 
spectra of higher shells to be fiatter than lower shells at the same rescaled 
frequencies. We thus conclude that dynamic scaling is incompatible with 
our numerical results. In the next section we will analyze the GOY model 
analytically and show explicitly where dynamic scaling fails. 

We also measured the cross correlations between shells. Recall that the 
symmetry (3) implies a vanishing of cross correlations that have a non zero 
modulo three difference between them. This effect was clearly observed 
in our simulation. Cross correlations of shells differing by a multiple of 
three were larger by two orders of magnitude compared to cross correlations 
of nearby shells without this property. The normalized cross correlations 
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decay as \n — m\ increases (by multiples of three) 
but seem to converge to a non zero value when one of the shells is in the 
dissipation range (larger than 33). 

Note that the period three oscillations observed by many authors [3] [6] 
are a direct result of aliasing of the effect of the forcing into the inertial range 
by the period three level to level correlations. They should be interpreted 
as a peculiar artifact of the GOY model and have no relation to 'lacunarity' 
in real turbulence. This is clearly shown by the fact that the period three 
behavior is independent of the value of the shell separation and will persist 
even in the limit of shell separation going to zero [7] [8]. On the other hand, 
the period 2 oscillations in the structure functions near the top of the inertial, 
are genuine dynamical effects, and are due to the fact that L is not positive 
definite [3]. 

In addition we measured the scaling of the structure functions 
These quantities are known to exhibit muUiscaling [2] and we have recovered 
this behavior. In fact, as we are using shells separated by \/2, we found 
stronger effects of intermittency (similar to what was observed in [6] ) . 

4 Direct interaction approximation for the 
GOY model 

As we have seen in the previous section, dynamic scaling fails to explain 
the shape of the frequency power spectrum. In order to examine this effect 
analytically we performed a resummed perturbation expansion for the GOY 
model. Basically we have derived the analog of the DIA equations [9] for the 
GOY model. The DIA approximation is known to fail for the Navier-Stokes 
equations due to the failure to respect Galilean invariance, or not taking into 
account the effect of large scale advection in a proper way [10]. The GOY 
model forms a truncated set of equations which are only coupled locally, 
hence conventional wisdom (which we will have no reason to doubt in what 
follows) indicates that the DIA should be a good approximation for the GOY 
model. 

As we have noted in section 2.2 the phase symmetries (3) induce the 
vanishing of some cross terms in the correlation functions and in the Green 
functions, however, not all of the cross terms vanish. This implies that, 
as opposed to the Navier-Stokes case, the GOY model DIA equations are 
formulated in terms of functions that are not diagonal in the shell n index, 
namely 

G{n,m,uj)6{uj + u') = {6un{oj)/6f^{uj')), 

(15) 

The non-diagonality makes for a slightly more complicated form for the 
DIA equations. They now contain expressions for off diagonal quantities as 
well as diagonal ones. It also complicates the discussion of the convergence 
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properties of the integrals and the sums. Nevertheless, one obtains the self 
consistent Dyson- Wyld integral equations, 



G{n,m,Lo) = G^{n,Lo){8n,m + ^^{n,l;to)G{l,m;Lo)) (16) 
U{n,m-Lo) = J2G*{n,k;u){{rji) + ^k,l;u))G{l,m;u), (17) 



k,i 

and the DIA for the mass operators 



i^m — j, —io — id') 



$(2)( 



n;(jj) = j d(jj'knk.raY^Y2^'''^33"i ^'^("^ ~^ 

ii' jj' 

xG{n + i',m- j',uj') (18) 
n;id) = ^ J duj'knkm c,i'Cjj'U{n + z, m + j, -u - u') 

ii' jj' 

x[/(n + z',m + j',cu'), (19) 
where the coupling coefficients are 

Cl,2 = C2,l = a, (^1,-1 = = fe, C_2,-l = C_i^_2 = c. 

Although in principle one should solve these equations to find the behavior of 
the correlation function, one usually tries a scaling ansatz in order to extract 
information about the scaling behavior of the correlation function U . Not 
withstanding our knowledge based on numerics that dynamic scaling fails, it 
is instructive to try a dynamic scaling form for G and U in order to see why 
it fails. We assume 

G{n,m,L0) = {knk.ra)~'^^g{^,k.ra/kn). ^^^^ 

U{n,m,uj) = {knkm)y/^f{j^,km/kn). 

where y and z are scaling exponents to be self consistently determined. If we 
assume naively that all sums and integrals converge in the limit of infinites- 
imal dissipation[ll] [p — )• 0) we find that the DIA equations for F and G 
both yield the same exponent relation 

z-y = 2 (21) 

Thus up to this point dynamic scaling gives us a one parameter family of 
solutions that is self consistent. In order to obtain a second exponent relation 
it is usual to use the conservation of energy in the inviscid limit (Euler equa- 
tion) to deduce a conservation of fiux in the inertial range (this conservation 
may be considered as the definition of the inertial range). The energy fiux 
through shell n, is defined in the GOY model as [6] [3] 

Cn = -'^m{-ckn-lUn-lUnUn+i + aknUnUn+iUn+2) ■ (22) 
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Assuming dynamic scaling, the lowest order terms in the expansion for the 
3-point correlators in (22) scale as kl^^'^'^^ ^ thus constancy of the the flux 
(tested numerically by [6]) implies that 

z + 2y = -2, (23) 

which with (2f) gives z = 2/3, y = —4/3, which is what a straightforward 
application of Kolmogorov type arguments gives [f]. 

Renormalized perturbation expansions in which Galilean invariance is 
preserved [f2] [f3] predict Kolmogorov scaling for Navier-Stokes turbulence, 
including dynamic scaling for the correlation functions and Green's functions. 
However, experiments indicate that the scaling of higher order structure func- 
tions differs from the K4f predictions. How to understand these deviations 
and how to reconcile these with the results of renormalized perturbation the- 
ories are a subject of active study [f4][f5]. At any rate, none of the present 
explanations resorts to assuming the failure of dynamic scaling as a cause. 
Indeed, in Navier-Stokes turbulence there seems to be no reason to doubt 
the existence of dynamic scaling. 

The GOY model is different in this respect. In addition to relationship 
(22) we have an additional conserved quadratic quantity L defined in 2.2 . 
This induces an additional conserved fiux 



in = -'<5m{{a/c){-kn-lUn-lUnUn+l + knUnUn+lUn+2)) ■ (24) 

which yields the exponent relation 

log(a/c) 

z + 2y = -2- ' ^ 25 
log 5 

in the same manner as (23) was obtained. The RHS of this scaling relation is 
not real in general. Such a solution is unacceptable, since the simultaneous 
pair correlation function must be positive. A more general scaling ansatz 
that includes a periodic function of n -|- m can solve this problem, and then 
only the real part of eq. (25) should be taken. However this scaling 

log \a/c\ 

z + 2y = 2+ f ' ' 26 
log 5 

contradicts the scaling relation (23) derived from energy fiux conservation. 
Indeed, the energy fiux conservation is violated when the original scaling 
ansatz is multiplied by a non-constant periodic function of n -|- m. 

There are two ways to solve the contradiction. Either one of the fiuxes 
vanishes in the inertial range (similar to the case of 2 dimensional turbulence) 
or we must reject the assumption of dynamic scaling. If one of the fiuxes 
is small in some sense compared to the other one may treat it as yielding a 
correction to dynamic scaling superimposed on top of the larger fiux's con- 
tribution. Our numerics do not support this case but we cannot rule out 



8 



the possibility that intermittency disappears asymptotically in the number 
of shells. The existence of two non vanishing fluxes in the inertial range is 
incompatible with dynamic scaling. A numerical check of the flux £ indeed 
demonstrates that it doesn't vanish. When only one level is forced, the two 
fluxes are linearly dependent [3], and therefore must be non-zero simultane- 
ously. We identify this fact as the reason for the failure of dynamic scaling 
observed numerically in section 3. 

5 A shell model with a single invariant 

Following the reasoning of the previous section, modifying the GOY model 
so as to remove the spurious L invariant should restore the possibility of a 
dynamic scaling solution for the DIA equations. A dynamic scaling solu- 
tion of the DIA entails the Kolmogorov scaling for simultaneous two point 
correlations. If the DIA is a good approximation then we cannot find a dy- 
namic scaling solution of the GOY with intermittency. This prediction can 
be tested in a model that obeys dynamic scaling, thus we must construct a 
model with only one conserved fiux. 

A simple way of eliminating the L invariant is by coupling the shell vari- 
ables to higher and lower shells. A generalization of the GOY interaction 
term that couples to more modes reads 

^ ( ^ ^ ^ij^n^n-\-i^n-\-j ~\~ ^ij^n — i^n — i^n-\-j — i ~\~ ^ij^n—j^n—j^n-\-i—j) • ('^'^) 
i<j 

The sum runs on a finite set of pairs (the original GOY is the special case 
where there is only one pair, namely {I, 2}). Energy conservation in the 
inviscid unforced case demands that 

«u + % + '^u = 0- (28) 

The condition (8) whose solutions are the quadratic invariants of the system 
is generalized in this case to 

Antt.j + An+tkj + An+jC,j = 0, for cvcry pair {ij}. (29) 

The condition (28) ensures that (29) is satisfied by An = const. The condi- 
tion (29) is clearly more restrictive than the original condition (8), and we 
will show in the following an explicit example in which energy is the only 
conserved quadratic invariant. 

In principle the pairs {ij} may chosen arbitrarily. However, if we want to 
ensure the vanishing of directly from the symmetry of the equations, we 
need the phase symmetry (3). Thus we are interested in including only pairs 
which preserve (3). This set includes all the pairs in which i -\- j is divisible 
by 3, and neither i nor j is divisible by 3. 

We now have quite a large parameter space at our disposal, of which we 
chose the next most local allowable pair which is {2, 4} in addition to {f , 2}. 
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Returning to equation (29) with the specific choice of pairs {2, 4}, {1,2} 
we see that there are two independent ways to obtain An+4 from 
One way, by iterating the equation for the {1,2} pair twice, and another 
way is by iterating it once to get an expression for An+2 and inserting this 
into the equation for the {2,4} pair. As both expressions must be equal, 
we find a mapping that gives An function of A„. As we know that 

An = An+i = An+2 = • • • is a solutiou due to (28), we conclude that there 
are no more independent quadratic invariants for this coupling choice. 

Our choice of couplings means that there is a 3 dimensional parame- 
ter space, since each pair involves 2 independent coupling constants and in 
addition a global scale factor can be absorbed in a redefinition of the Un 
variables. As observed by [16], the original GOY doesn't exhibit chaotic be- 
haviour in parts of the parameter space; instead there can be a stable fixed 
point, limit cycle, or quasi-periodic motion. This phenomenon persists also 
in the enlarged parameter space [8]. We expect dynamic scaling to appear for 
parameter values for which there are chaotic solutions. After a qualitative 
mapping of the parameter space we found a region which has the desired 
property. Typical parameter values in this region are 

ai2 = 1 bi2 = -0.5 ci2 = -0.5 ^^^^ 

(224 — 5 ^24 — — 2.5 C24 — — 2.5 

This choice of parameters is essentially two original GOY-like couplings, the 
second of which is 5 times stronger than the first. The fact that the {2,4} 
coupling is much stronger than the {1,2} coupling raises a problem that if 
only one level is forced, the shells with the same parity as the forced shell 
will have larger amplitudes than neighboring shells with opposite parity. This 
problem is dealt with by forcing shell 2 by /2 = (1 -|-z)10~^ in addition to the 
usual forcing of shell 3, but even then there remain some period 2 oscillations 
in the structure functions. Other numerical parameters remain unchanged, 
including the simulation duration. 

We have measured the same quantities as for the original GOY model 
(3). The data collapse for the frequency domain power spectra is shown in 
figure 2. The fit is clearly superior to the fit accomplished in figure 1. Also, 
the empirical rescaling parameters hn and cOn defined in 14 can be fitted well 
to power laws 

COn ~ kf^-'' (31) 
hn ~ k-'-''^-'^ (32) 

Which are compatible with the dynamic scaling prediction and are close to 
the Kolmogorov exponents 2/3,-4/3. As another check we have measured 
higher order structure functions and compared their scaling exponents 

to those of the original GOY and to Kolmogorov scaling. Table 1 contains 
this comparison. Note that higher moments still do not conform to the 
Kolmogorov exponents but the deviation is much smaller. 
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The deviation may be attributed to many sources, such as finite system 
efi'ects (corrections to scahng), a persistence of some memory of the extra 
invariant and, of course, real intermittency, indicating a failure of the per- 
turbation scheme on the level of the DIA. It is beyond the scope of this paper 
to decide which of the above effects is responsible for the deviations from Kol- 
mogorov scaling. Our purpose is to stress that these effects are much smaller 
than the ones observed in the original GOY model. We attribute the large 
deviations in the original GOY to the lack of dynamic scaling that is induced 
by the existence of a second fiux in the inertial range. 

6 Conclusions 

Numerical evidence shows that the solutions of the GOY model do not have 
dynamic scaling. In addition, the DIA equations do not have solutions which 
exhibit dynamic scaling, due to the existence of two non vanishing constant 
fiuxes in the inertial range. The derivation of Kolmogorov scaling for the 
structure functions via perturbation methods relies on the existence of dy- 
namic scaling. Therefore, there is no evidence for a contradiction between 
perturbation theory and multiscaling in the GOY model. 

This situation is unlike the case of Navier-Stokes turbulence. Experimen- 
tally, high order structure functions deviate from the predicted Kolmogorov 
exponents [17]. However, we have no good reason to doubt the scale in- 
variance in time and space of turbulent motion. Helicity, or other conserved 
quantities which in principle give rise to conserved fiuxes, are zero in isotropic 
homogeneous turbulence, and are considered small in theories of turbulence. 
A dynamic scaling ansatz in the perturbation theory is then consistent and 
yields the Kolmogorov scaling uniquely [12] [18]. It should be stressed that 
dynamic scaling has never been proved from first principles, and its failure 
may be an interesting alternative mechanism for explaining some aspects of 
turbulent phenomenology. The results of this paper motivate us to propose 
that dynamic scaling should be tested directly in experiments in turbulence. 
A major deviation from dynamic scaling would imply that a reconsidera- 
tion of the current working assumptions in theories of turbulence is needed. 
Whether or not the final resolution of intermittency in turbulence will rely 
on presently applied techniques or whether non-perturbative effects should 
be considered is still not clear. 

The reason for the difference between the Navier-Stokes intermittency and 
the intermittency in the GOY model is the existence of a second quadratic 
invariant. Intermittency in a model in which the second invariant is absent 
may have an origin closer to the Navier-Stokes type intermittency. However, 
our results seem to indicate that a modified GOY model without the second 
invariant has only small remnants of intermittency. Therefore, in our opinion, 
the relevance of multiscaling in the GOY to the dynamics of Navier-Stokes 
is questionable. 
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n 






1 


2 


3 


4 


5 


6 


Original GOY 


^q= 




.42 


.78 


1.05 


1.20 


1.29 


1.35 


Modified GOY 


(q= 


=V2) 


.36 


.69 


1.00 


1.28 


1.53 


1.81 


Original GOY 


^q= 


=2) 


.39 


.75 


1.07 


1.36 


1.64 


1.87 


Modified GOY 


(q= 


=2) 


.35 


.69 


1.03 


1.35 


1.65 


1.92 


K41 






.33 


.67 


1.00 


1.33 


1.67 


2.00 



Table 1: A comparison between the scaling exponents (n of the structure 
functions in a few models and Kolmogorv scaling. The errors in the scaling 
exponents are approximately ±.01, except for the modified model with q = 2 
where they are ±.03, due to large fiuctuations in the structure functions. 




Figure 1: Frequency power spectra for three shells 12(dashed), 22(full), 
32(dot dashed) after the rescaling given in text. The intersection point is 
an artifact of the rescaling and not a real feature of the data. Although there 
is a superficial similarity, the three functions depicted in this figure differ sig- 
nificantly. We measured the difference by calculating the ratios (c<;")/(c<;^)"/^ . 
These ratios deviate significantly between shells indicating that the functions 
do not collapse onto a universal scaling function. 



13 



10.0000 E 




Figure 2: Same as figure 1 for tlie modified amodel. Note tlie significant 
improvement in tlie fit between tlie rescaled power spectra. 
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